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Abstract 

Using an innovative technique arising from the theory of symmetric spaces, we obtain an 
approximate analytic solution of the Dorokhov-Mello-Pereyra-Kumar (DMPK) equation 
in the insulating regime of a metallic carbon nanotube with symplectic symmetry and an 
odd number of conducting channels. This symmetry class is characterized by the presence 
of a perfectly conducting channel in the limit of infinite length of the nanotube. The deriva- 
tion of the DMPK equation for this system has recently been performed by Takane, who 
also obtained the average conductance both analytically and numerically. Using the Jaco- 
bian corresponding to the transformation to radial coordinates and the parameterization 
of the transfer matrix given by Takane, we identify the ensemble of transfer matrices as 
the symmetric space of negative curvature SO*(Am + 2)/[SU(2m + 1) x £7(1)] belonging to 
the Dill-odd Cartan class. We rederive the leading-order correction to the conductance 
of the perfectly conducting channel (In 8g) and its variance Var(ln 5g). Our results are in 
complete agreement with Takane's. In addition, our approach based on the mapping to 
a symmetric space enables us to obtain new universal quantities: a universal group the- 
oretical expression for the ratio Var(ln5g)/ (ln5g) and as a byproduct, a novel expression 
for the localization length for the most general case of a symmetric space with BC m root 
system, in which all three types of roots are present. 



1 Introduction 



Carbon nanotubes, consisting of single or multiple graphene layers wrapped into nanometer- 
thick cylinders 1 , have been an object of experimental and theoretical study in mesoscopic 
physics over the past decade. They exhibit particular transport properties due to the spe- 
cial characteristics that distinguish them from ordinary quasi-one dimensional quantum 
wires. The hexagonal lattice of carbon atoms puts constraints on the possibilities of match- 
ing a single layer of atoms onto itself to form a cylinder or torus, giving rise to systems with 
differing discrete chirality, and in the case of a toroidal geometry, differing degree of twist. 
The cylindrical or toroidal geometry allows the conduction electrons to orbit around the 
axis of the cylinder. This happens for certain chiralities only. Experimental and theoretical 
analyses show that carbon nanotubes are either metallic or semiconducting. The electrical 
properties of carbon nanotubes are governed by electrons close to the Fermi energy and 
show some interesting phenomenology (see e.g. [T]). 

A carbon nanotube can be considered a special kind of disordered quasi-one dimensional 
quantum wire, and the same techniques that are applied in the theoretical study of meso- 
scopic wires can be applied also to carbon nanotubes. Disordered systems were classified 
by Dyson into three standard symmetry classes (orthogonal, unitary and symplectic, la- 
belled by Dyson index (3 = 1, 2, and 4 respectively) 0. In this paper we will focus 
on metallic carbon nanotubes with symplectic symmetry. By this we mean that the sys- 
tem is time-reversal invariant but has strong spin-orbit scattering, according to Dyson's 
classification. 

Such nanowires have been studied by Ando and Suzuura who showed that they have 
an anomalous transport property: there is one perfectly conducting channel present even 
in the limit of an infinitely long wire. The first discovery of the non-vanishing conduc- 
tance in the long-wire limit in the symplectic symmetry class is, however, due to Zirnbauer 
[3] and to Mirlin, Miiller-Groeling, and Zirnbauer jH]. These authors expressed moments 
of the conductance of thick disordered wires in terms of the heat kernel in a supersym- 
metric a-model approach, subsequently obtaining approximate analytical expressions for 
the conductance for all three symmetry classes by Fourier analysis on the supersymmetric 
manifold. 

Takane jEl El El and Sakai and Takane $\ conducted further studies of systems with sym- 
plectic symmetry and with an even or odd number of conducting channels, employing a 
variety of methods. It was confirmed in jB] and [7j, using the supersymmetric formalism 
and numerical simulations, respectively, that for an odd number of conducting channels 
(which is the case in a carbon nanotube), the dimensionless conductance g — > 1 in the 

1 Stacked sheets of single graphene layers form graphite. 
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limit of an infinitely long wire, indicating the absence of Anderson localization. In this 
case there is an unpaired transmission eigenvalue equal to unity which gives rise to the 
anomalous behavior of the conductance. This is to be contrasted with the case of an even 
number of conducting channels, when the average dimensionless conductance approaches 
zero with increasing wire length L (this case is realized in an ordinary quantum wire, where 
the limit L — > oo corresponds to the insulating regime). 

In |H| these results were confirmed using scaling theory. Within the framework of random 
matrix theory, the Dorokhov-Mello-Pereyra-Kumar (DMPK) equation for the probability 
distribution of the transmission eigenvalues as a function of the length of the wire was 
constructed for this particular system, and the conductance and its variance were derived 
in the limit of an infinitely long wire. In contrast, it was proved in [9] that the presence 
of the perfectly conducting channel in the odd-channel case does not have any bearing on 
the conductance and its variance if the length of the wire is shorter than the localization 
length, in which case there is almost no even-odd difference. This is due to geometric 
eigenvalue repulsion in random matrix theory. As a consequence of the presence of the un- 
paired transmission eigenvalue 2 T 2m+ i = 1, the other eigenvalues are repelled and therefore 
diminish (we recall that < T 4 < 1). 

This paper is organized as follows. Sections El and El contain some generalities concerning 
the symmetric space description of random matrix theories and the standard theoretical 
description of mesoscopic scattering, respectively. In section |U we discuss the mapping of 
the DMPK equation to the equation of free diffusion on the symmetric space. 

In section |S] we then identify the symmetric space corresponding to the ensemble of transfer 
matrices for the symplectic ensemble with an odd number of scattering channels. This will 
be done in two ways: first using the Jacobian on the symmetric space, and then using the 
parametrization of the transfer matrix given by Takane jH] . In section |B] we show how to 
solve the DMPK equation and obtain the conductance of the system using the technique 
of zonal spherical functions on the symmetric space. Some details of the calculation are 
reported in Appendix |XJ 



2 The symmetric space point of view 

The classification of random matrix ensembles in terms of symmetric coset spaces of Lie 
groups (or as subspaces of Lie algebras) is by now rather well-known. In reference jTUj 

2 We use Takane's notation in which N = 2m in the even-channel case and N = 2m + 1 in the odd- 
channel case. In this notation, the number of non-zero, distinct, paired eigenvalues in the symplectic 
symmetry class is m. (This quantity is often denoted by N elsewhere.) 
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a detailed review of this subject was given. It was shown that all the known hermitean 
random matrix ensembles can be identified with symmetric coset spaces of Lie groups or 
with subspaces of Lie algebras. There are exactly twelve 3 such groups of symmetric spaces. 
Each group comprises three symmetric spaces identified with manifolds of positive, zero, 
and negative curvature, respectively. The positive curvature spaces are identified with 
ensembles of scattering matrices, the zero curvature spaces (which are Lie algebra sub- 
spaces) with ensembles of Hamiltonians, and the negative curvature spaces with ensembles 
of transfer matrices. 

In the third and fourth column of Table ^ which we have reproduced from reference ^0] , 
the symmetric coset spaces of positive and negative curvature arising from simple Lie 
groups are listed (the zero curvature spaces are not listed since these are simply identified 
with the Lie algebra subspace defining the corresponding curved manifolds). In the fifth 
through seventh column of Table Q the multiplicities of the roots of the corresponding 
restricted root systems are listed. The subscripts refer to ordinary, long and short roots, 
respectively. 

Let us remind the reader of the construction of the root system of a simple Lie algebra. 
There are two kinds of generators in the algebra: there is a maximal abelian subalgebra 
{Hi, H r }, where r is the rank of the algebra and [Hi,Hj] = 0, and there are raising 
and lowering operators E± a . The roots a are functionals on the Cartan subalgebra whose 
components «j satisfy 

[H i ,E a } = a i E a (2.1) 

The angle between root vectors can take only a few discrete values, which makes Cartan's 
classification possible and leads, for the non-exceptional Lie groups, to four infinite series 
of root systems A„_i, B n , C n , and D n . Just as a root system corresponds to each simple 
Lie group, for each symmetric coset space one can identify a unique restricted root system 
where each root vector has a given multiplicity. It would bring us too far to discuss the 
construction of these root systems here, and we refer the interested reader to JU] for details. 
In any case it is similar to the construction of the root system for the entire Lie algebra. 

The non-reduced root system BC n is defined as the union of B n and C n . It is of relevance 
in this context as it appears as the restricted root system of several types of symmetric 
spaces, as can be seen by inspection of Table [T] 

3 In Cartan's classification there were eleven but it is natural to split the DJ/I-symmetry class into 
even and odd, following the example of some recent authors. This is because the even and odd cases 
correspond to different restricted root systems. 
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Table 1: Cartan's classification of irreducible symmetric spaces and some of their random 
matrix theory realizations. In the third and fourth column we have listed the symmetric 
coset spaces of positive and negative curvature based on simple Lie groups. The zero cur- 
vature spaces are not listed since these are simply identified with the Lie algebra subspace 
defining the corresponding curved manifolds. m Q , mi, and m s are the multiplicities for 
the respective restricted roots and in the last three columns we find the random matrix 
ensembles with known physical applications corresponding to symmetric spaces of positive, 
zero and negative curvature, respectively. The ensemble of transfer matrices for symplec- 
tic odd-channel carbon nanotubes corresponding to the present case has been inserted in 
boldface in the last column. We have set v = p — q and for the random matrix ensem- 
bles we use the following abbreviations: C for circular, G for gaussian, \ f° r chiral, B for 
Bogoliubov-de Gennes, P for p-wave, T for transfer matrix and S for S-matrix ensembles. 
The upper indices indicate the curvature of the respective symmetric spaces, while the 
lower indices correspond to the multiplicities of the restricted roots characterizing each 
triplet. 
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In this framework, the multiplicities of restricted root vectors associated to a symmetric 
coset space are identified in the random matrix ensemble as the exponents appearing in the 
Jacobian for the transformation from the space of random matrices to the space of random 
matrix eigenvalues. This Jacobian gives rise to the characteristic eigenvalue repulsion in 
random matrix theory. Such a transformation may be expressed as 

S -> U^XU (2.2) 

where S is a random matrix, X is (block-) diagonal and contains the random matrix eigen- 
values and U is unitary. The random matrix partition function is invariant under such 
a transformation. For example, for the simple gaussian unitary random matrix ensemble 
characterized by root multiplicities m a = 2, mi = m s = 0, the Jacobian of such a coordi- 
nate transformation is given by J({xi}) ~ Y\i<j( x i ~ x j) 2 (where {xi} are the eigenvalues). 
In this case there is only one exponent in the Jacobian, (3 = 2 (where (3 indeed is Dyson's 
symmetry index) which is given by the multiplicity of ordinary roots: (3 = m Q = 2. Sim- 
ilar relations between the exponents in the Jacobian and the three root multiplicities are 
true also in the case of more complicated Jacobians with more than one exponent. We 
conclude that the root structure of the underlying symmetric space completely determines 
the celebrated universal random matrix eigenvalue correlations! 

The random matrix ensembles identified with the symmetric spaces in each row are listed 
in the last three columns in Table H The notation for the various types of ensemble is 
as follows: C for circular, G for gaussian, \ for chiral, B for Bogoliubov-de Gennes, P 
for p-wave, T for transfer matrix and S for S-matrix ensembles. A discussion of these 
ensembles here is outside the scope of this paper and we refer to ^U] or to the vast random 
matrix literature for an introduction. 



3 Scattering in a mesoscopic wire 

All the material in this section is standard and is included only to make the paper self- 
contained. 

The scattering of electrons in a mesoscopic wire attached to ideal leads can be described 
using the transfer matrix. Assuming there are m propagating modes at the Fermi level, 
we describe them by a vector of length 2m of incoming modes /, I' and a similar vector 
of outgoing modes O, O' in each lead. Let the unprimed letters denote the modes in the 
left lead and the primed letters the modes in the right lead. While the scattering matrix 
S relates the incoming wave amplitudes /, I' to the outgoing wave amplitudes O, O', the 
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transfer matrix M relates the wave amplitudes in the left lead to those in the right lead: 





(3.1) 



The transmission eigenvalues are the m eigenvalues of the matrices tv (or equivalently, 
ft''), where t (f) is the left-to-right (right-to-left) transmission matrix appearing in the 
scattering matrix (and r, r' are the respective reflection matrices): 



S ={tr>) < 3 - 2 > 
The transfer matrix is expressed in terms of these submatrices as ^T] 

M-(%KV) (3.3) 



In terms of the transmission eigenvalues {Tj}, the expression for the dimensionless conduc- 
tance is given in the Landauer-Lee-Fisher theory by 



G m 2e 2 

Cr j =1 II 

Instead of the transmission eigenvalues, it is common to use the non-negative variables 
{Xi} defined by 



A, = ^ (3.5) 

M can then be parametrized as 4 |12j 

4 This parametrization is valid for the cases M E Sp(2m ) R)/[SU(m) x U(l)] ("orthogonal" ensem- 
ble) and M £ SO* (4m)/ [SU (2m) x U(l)] ( "symplectic" ensemble), whereas for the unitary ensemble 
SU(m, m)/[SU(m) x SU(m) x 1/(1)], u* and v* have to be substituted by two unitary matrices v! and v' 
not related to u and v by complex conjugation. 
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where u, v are unitary m x m matrices and 

A = diag(Ai, A m ) (3.7) 

is a diagonal matrix. 

In the symplectic (even-channel) case, there is a doubling of the degrees of freedom due 
to the fact that the components of I, O, etc. become spinors. This means the eigenvalues 
{A,}, i = 1, ...,m are doubled so that we get m pairs of degenerate eigenvalues (this is 
referred to as Kramers degeneracy). If we write u, v and A in terms of complex matrices, 
we are dealing with 2m x 2m matrices, whereas in terms of quaternion matrices, they are 
m x m. The number of distinct eigenvalues is m. In this case eq. ()3.6|) and ()3.7|) are still 
valid, but the matrix elements are now quaternions. Following Takane, in the even-channel 
case we then set A" = 2m and in the odd-channel case A" = 2m + 1, so that the matrix M 
in the symplectic ensemble becomes a matrix with N x N complex elements. 

It is known that the parametrization ()3.6|) leads to a coset space structure for A" = 2m 
|13| I14j. One can show that the symmetries - in this case flux conservation, time reversal 
symmetry, and no spin-rotation invariance - imposed on M lead to M G SO* (4m) for 
the symplectic even-channel case, in which all the eigenvalues are paired due to Kramers 
degeneracy and there is no zero eigenvalue. It is easy to see that the eigenvalues {Aj} are 
unchanged if one transforms the matrix M in the following way: 

M -> M' = WUTVW^ 1 = U'YV (3.8) 

where W for the symplectic ensemble is a unitary quaternion matrix. Thus, in the symplec- 
tic even-channel case, M belongs to the symmetric coset space SO* (4m) /[SU (2m) x £7(1)] 
corresponding to the Dill-even Cartan class. 
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4 DMPK equation as free diffusion on a symmetric 
space 



The DMPK equation is a differential equation describing the evolution of the transmission 
eigenvalues of a quantum wire with increasing length of the wire. Its solution gives access 
to the distribution of transmission eigenvalues and therefore to the conductance through 
equation (JH3J). 

In the random matrix approach to quantum wires, as a consequence of the description of 
random (transfer) matrix ensembles in terms of symmetric spaces, the DMPK equation 
of a quantum wire can be expressed in terms of the radial part of the Laplace-Beltrami 
operator on the appropriate symmetric space. The Laplace-Beltrami operator (denoted 
Ab) is simply the operator of free diffusion on a given symmetric space, and corresponds 
to the Laplacian in R 3 . It is defined as the lowest order Casimir operator of the algebra, 
expressed in local coordinates as a differential operator. We denote its radial part by 
A^. As an example, the Laplace-Beltrami operator on the symmetric space SO(3)/ SO(2) 
(whose points are in one-to-one correspondence with the two-sphere) is given, in radial 
coordinates, by Ab = <9J + cot# dg + sin -2 ^ (we have set the radius of the sphere r = 1), 
and its radial part is A' B = dg + cotOdg. 

It is a completely general result, discussed for example in [T3J [TU], that the DMPK equation 
for the transmission eigenvalues can be written as 

dP({ ^ hs) =(2 1 )- 1 BP({x t },s) (4.1) 
where 7 is a constant and {xi} are radial coordinates on the symmetric space. The operator 



B = t4rJ(M)^-J-\M) (4.2) 

k=l L/- t/ fc 

is related to the radial part of the Laplace-Beltrami operator A' B on the symmetric space 
through 



A' B = r l BJ (4.3) 
and J({xj}) is the Jacobian for the transformation from the space of random matrices to 
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the space of eigenvalues (radial coordinates), cf. section |21 The constant 7 is in the most 
general case given by (cf. the expression in reference [14J where 9 = 0) 



2(1 + ri + /3(m-l) + 9/2) 
7 I + 77 ' 

where 77 = mi, (3 = m Q , 9 = m s and mi, m a , m s denote the multiplicities of the long, 
ordinary, and short (restricted) roots of the symmetric space and f3 is Dyson's index. This 
expression was obtained simply by matching 7 to the explicitly known cases. In eq. ()4.1|) 
the dimensionless variable s is defined by s = L/l, where L is the length of the mesoscopic 
wire and / is the mean free path of the scattered electrons. 



5 Determination of the symmetric space 



In this section we will determine the symmetric space to which the transfer matrix of the 
symplectic odd-channel carbon nanotube belongs, using two different methods. 



5.1 Determination using the Jacobian 

The DMPK equation for the symplectic ensemble with an odd number iV = 2m + 1 of 
conducting channels was derived in jH] and found to be given by 

8P(\ u ...,\ m ,s) 1 " 9 / ;8%^i)\ 

* = /V^S9A-( A - (1 + A ' )J 9A- J J (5 ' 1) 

where we believe there is a typing mistake in a sign in the formula given in jH]. Note that 
there are m pairs of distinct non-zero eigenvalues and one zero eigenvalue in this case. For 
an odd number of conducting channels, J is given by 



m m 



^({A l }) = n^n(^-A fc ) 4 (5.2) 

i=l j>k 

In jH], the corresponding function was given for the even-channel case: 
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m 

£({*}) = IK**-**) 4 

j>k 



(5.3) 



The coordinates {Aj} are not the radial coordinates on the symmetric space. By performing 
the variable substitution 

Aj = sinh 2 Xj (5.4) 

in the integration measure 5 

km) n dA * Ufa}) n dx i ( 5 - 5 ) 

i i 

and using identitites for the hyperbolic functions, we obtain the Jacobian on the symmetric 
space, expressed in the radial coordinates {xi}\ 

mm m 

Jo({xi}) = Y\_ sinh 4 Xj Y\_ sinh 4 (xj — Xfc)sinh 4 (xj + Xt) Y\_ sinh(2xz) (5.6) 

1=1 j>k 1=1 

Using now the general result for the Jacobian of the transformation to radial coordinates 
on a symmetric space, discussed in detail in [TU] . 

J ( ' ] (x) =n aeR+ (sinh(x a ))^ (5.7) 

where 0, — , + denote the curvature of the space, R + denotes the set of positive restricted 
roots, x a denotes the projection of x on the root a, and m a its multiplicity, we conclude 
that eq. (J5.6)) describes the Jacobian on a symmetric space with negative curvature and 

5 It is important to note that J ({A;}) \\ i d\ is the integration measure on the symmetric space ex- 
pressed in the "wrong" coordinates. Simply doing the substitution (|5.4|l in J„({Ai}), we would miss the 
factor coming from the long roots. This is evident in going from eq. I|4.1|l to eq. I|5.1|l . where we have to 
take into account also the factors of -g- 2 - in the chain rule expression for . 
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root multiplicities m s = 4, m a = 4, mi = 1 for the short, ordinary, and long roots, 
respectively. Comparing with Tabled we conclude that the transfer matrix M in the odd- 
channel case belongs to the symmetric space SO*(4m + 2)/[SU(2m + 1) x U(l)] (Cartan 
class Dill-odd). We have inserted this ensemble in the last column of Tabled in boldface. 

A similar analysis of the Jacobian for the even-channel case using eq. (|5.3jl results in the 
coset space SO* (4m) /[SU (2m) x U(l)] (D Ill-even), which agrees with the known result 
from the parametrization ()3.6|) given above. 



5.2 Determination using a parametrization of M 

The symmetric coset space of M in the odd-channel case can also be identified using the 
parametrization 



given by Takane jHJ. Here 9 = ih where h is an (2m + 1) x (2m + 1) hermitean matrix and 
r\ is an (2m + 1) x (2m + 1) arbitrary complex antisymmetric matrix. 

Let us compare this with the even-channel case. In ^T] it was observed that in the even- 
channel case, the parametrization (J3.6)) is equivalent to 



with 6 defined as above except that it is now a 2m x 2m matrix, and £ is an arbitrary 
2m x 2m complex symmetric matrix, if one makes the identifications 




(5.8) 




(5.9) 



e = uv ( = v 



(5.10) 



Note that this agrees with the properties of 9 and £, if \/~K. is symmetric: 
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\ 



A m / 



(5.11) 



where each eigenvalue is twofold degenerate. It is easy to verify that Takane's parametriza- 
tion ()5.8|) in the odd-channel case is equivalent to ()3.6j) if instead of equation (|5.10j) we 
identify 



e e = uv i] = v^vA'v 



In,* 



(5.12) 



where u, v are unitary and v^/V denotes an antisymmetric (2m + 1) x (2m + 1) matrix of 
the form 



/ 0_ 

Ai 



V 



VA 

-VKt o 



(5.13) 



J 



Note that any antisymmetric matrix rj can be written in the form expressed in equa- 
tions (|5.12I5.13|) . Therefore we can use the same argument as in the even-channel case 
(eq. (j3.8j0 to derive the coset structure SO*(4m + 2)/[SU(2m+ 1) x U(l)] of the ensemble 
of transfer matrices in the odd-channel case. The restricted root system corresponding to 
this symmetric space is of the type BC m . This will be used in the following section, where 
we will also remind the reader of some general features of root systems. 

To our knowledge this is the first known physical realisation of this ensemble, which in the 
notation of table [T] we denote T± x 4 . Thus we have now filled out the previously empty 
space in Table ^ corresponding to this physical realisation. As was noted also in ^U] , the 
empty spaces in this table are not really "empty" , but so far no application of them known 
to us has been discussed in the literature. 
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6 Conductance from symmetric spaces 



Having identified the symmetric space to which the transfer matrix M belongs, we can 
now solve the DMPK equation using the technique of zonal spherical functions ^U]. In 
reference ^2] the exact procedure for solving the DMPK equation using the method of 
zonal spherical functions was described in some detail. Here we will briefly outline the 
steps, but the full details of the computation can be found in the Appendix. A similar 
computation in the exactly solvable (3 = 2 case can be found in [TE] . 

The zonal spherical functions are eigenfunctions of the radial part of the Laplace-Beltrami 
operator on a symmetric space: 

A> B Mx)=lA' B (k)Mx) (6-1) 

(where for brevity we have set x = {x±, x m }, k = {ki, k m } and 7a' (k) are known 
eigenvalues that may be functions of the roots, see for example [T7| or flUj ) . By using 
equation (J4.1j) and the mapping A' B = J~ 1 BJ where B is the DMPK operator given in 
eq. (|4.2j) . we see that if is a zonal spherical function corresponding to the radial part 

of the Laplace-Beltrami operator A' B on a given symmetric space, then J(x)<pk{x) is an 
eigenfunction of B/{2 r y) corresponding to the same symmetric space (and with the same 
eigenvalue multiplied by 1/(27)). 

The functions <fik{x) have been studied by Harish-Chandra [IB]. They are related to irre- 
ducible representation functions of groups and have a deep group theoretical significance 
[TO] ITS] . They form a complete basis in the space of square-integrable functions on the 
symmetric manifold, and can be used to express the analog of a Fourier transform on the 
symmetric space. Accordingly, the eigenvalue density we are seeking can be written in the 
form 

P(x, s) = J{x)f{x) = J(x) I f{k)<t> k {x)-^- 2 (6.2) 

where the integral over k defines a Fourier transform on the symmetric space of the function 
f(k) describing the initial conditions in fc-space for the solution P(x, s) of the DMPK 
equation. We choose f{ki) to be a gaussian: 

fih) oc e"^ (6.3) 
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where s and 7 have the same meaning as in the DMPK equation. This choice corresponds 
to ballistic initial conditions I ^> L [T^]. The expression for c(k) is given by 



c(k ) = n r(?(fc ' Q)/2) (6 4) 

where the product goes over the positive roots of the root lattice R, T is the Euler gamma 
function, m a the multiplicity of the root a, and (k, a) denotes a scalar product. 

The zonal spherical functions are known exactly for the ensembles of transfer matrices cor- 
responding to Dyson index (3 = 2. For Dyson indices /3 = 1, 4 only asymptotic expressions 
are known. We will use the expression for |c(fc)| 2 together with the asymptotically known 
large-x form of 4>h{x), valid for all values of k: 



4>k{x) ~ J- 1/2 (x) c(rk)e i{rk > x) (6.5) 

where rk is the vector obtained by acting on k with the element r of the Weyl group. We 
remind the reader that the Weyl group W is the symmetry group of the root system. It 
consists in the most general case of reflections and permutations of the root vectors. 

We are particularly interested in the solution to the DMPK equation in the insulating 
regime, where the peculiar feature of the carbon nanotube - the existence of one perfectly 
conducting channel in the long-wire limit - has the most relevant physical effect. In order 
to keep our analysis as general as possible, we will discuss the solution of the DMPK 
equation for a generic BC m root lattice (i.e. for a root lattice of rank m with general 
nonzero values for all three types of roots). As before, we use the following notation for 
the root multiplicities: mi = rj, m = (3 and m s = 9. 

Following the method outlined above, and described in full detail in Appendix^l we obtain 
for the probability density P({x«}, s) of the DMPK equation in the insulating regime 



P({x n }, s) oc ]J(sinh 2 Xj - sinh 2 Xi) 13/2 (x 2 - x 2 ) 

j>i 



x 

k 



where ^(77) = (h(9) = 0) if no long (short) root is present in the lattice (i.e. if 77 = 
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(9 = 0)), otherwise ^(77) = 1 (h{9) = 1). In the present case h{rj) = h{9) = 1. 



After ordering the x n 's from small to large and using the fact that in this regime 1 <C x± <C 
X2 <C ■ ■ • <C x m , we can approximate the eigenvalue distribution as follows: 



m 

P({x n },s)ocl[e-^-^ 2 (6.7) 

i=l 

where 



Xi=-(Tj + - + l3{i-l)) (6.8) 
7 2 

In the strongly localized limit the conductance is dominated by the first transmission 
eigenvalue. In this case, due to the presence of the perfectly conducting channel, the first 
eigenvalue actually gives the leading order correction to the contribution of this channel. 
We obtain the leading contribution Sg = from equation ()3.4|) . Using also ()3.5|) and 

(J5~4) we find: 

QnSg) = -y ~ InTi « -2x x = _^±ll ( 6 . 9 ) 
where we have defined the localization length £ which we find to be given by the expression 



i = -^j (6.10) 

The expression for 7 was given in equation (J4.4j) . Substituting the values of the root 
multiplicities in the present case (r) = 1, (3 = 4, 9 = 4) and defining (following [Uj) 
iV = 2m + 1 we find, in complete agreement with Takane's results, 



£=l(N-l)l (6.11) 

-a n ^)~(i) (6-12) 

Var(ln^) = ^ T (y) (6.13) 
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where Varx = (x 2 ) — (x) 2 is the variance of x. It is interesting to note that these two 
quantities can be combined into a universal ratio which does not depend on the microscopic 
details of the model (i.e., on the value of I or on the number of open channels) but only 
on group theoretical quantities: the root multiplicities of the appropriate symmetric space. 
Its general value is: 



Var ( ln ^) = 2 (6 14) 

(InSg) " r] + 6/2 1 ; 

This value only depends on the multiplicities of the short and long roots and could play 
the same role in the insulating regime as the universal conductance fluctuations play in 
the metallic regime (the latter depend only on /3, the multiplicity of the ordinary roots). 
In all the "standard" quantum wires this ratio is always 2, because in the corresponding 
ensembles the root multiplicities have the values rj = 1 and 6 = 0. This makes the present 
carbon nanotube case particularly interesting. Since in this case 6 = 4 it is the first 
example 6 of a nontrivial value of this universal ratio. 



7 Conclusion 



We have studied the electric conductance of single-layer carbon nanotubes with an odd 
number of conducting channels and belonging to the symplectic symmetry class, defined 
as the universality class of time-reversal invariant systems with strong spin-orbit scatter- 
ing. The conductance of metallic carbon nanotubes can be studied within the Landauer 
formalism used for ordinary quantum wires. To obtain the distribution of transmission 
eigenvalues needed to compute the conductance in this formalism, we have employed a 
technique related to symmetric spaces. 

After identifying the relevant symmetric coset space SO*(4m + 2)/[SU(2m + 1) x U(l)} 
to which the transfer matrix of the system belongs, we show how knowledge of the theory 
of symmetric spaces may be used to obtain a solution to the Dorokhov-Mello-Pereyra- 
Kumar (DMPK) equation, a differential equation describing the evolution of the trans- 
mission eigenvalues with the length of the carbon nanotube. The solution is obtained by 
mapping the DMPK equation onto the equation for free diffusion on the symmetric space 
defined by the ensemble of transfer matrices. We then find a suitable approximation to 
the distribution of eigenvalues in the (would-be) insulating regime, where the existence of 
a conducting channel even in the limit of an infinitely long nanotube in the odd-channel 

6 See [21 for another possible candidate, for which there is, however, not yet any clear-cut physical 
evidence. 
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case is known. Knowing the distribution of transmission eigenvalues, we can compute the 
leading contribution to the physical conductance of the carbon nanotube in terms of root 
multiplicities of the symmetric space. Our results can also be expressed in a form that is 
in complete agreement with Takane's results. 

In our case we have the maximal number of different types of roots (i.e. a nonzero number 
of long, ordinary and short roots in the BC m root system corresponding to the symmetric 
space manifold). As a consequence, we obtain the most general expression for the local- 
ization length of such a system in terms of root multiplicities (eq. (jfi.lOj) ). In addition we 
find the expression for the universal ratio 



in the insulating regime, which may play a role similar to that of universal conductance 
fluctuations. These universal group theoretical results are new and could only be derived 
using our technique based on symmetric spaces. 
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A Appendix 

In this Appendix we give all the details of the calculation of P({xi}, s). As we have already 
discussed, we can start from the expression involving the Fourier transform 



Var(ln 8g) 



2 



(7.15) 



(In Sg) 



V + 0/2 



P(x 7 s) = J(x)f(x) = J(x) I f(k)Mx) 



dk 



(A.l) 



c(k)\ 2 



where J(x) is the Jacobian on the symmetric space and 



f(ki) (x e~^7 



(A.2) 
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defines the initial condition. The function c(k) is given by 



c(k) - TT r(z(fc,q)/2) 

Explicitly, for a root lattice with long, ordinary and short roots with multiplicities mi = 

1, m = 4, m s = 4, c{k) is given by 



c(k) ft r( ' fci) ft r( ^ } ft F( ^ J (A 4) 

C{ ] \{ r(§ + }± r(2 + r(2 + i*±±) M r(2 + i%) 



We will be interested in the insulating regime in which the variables {ki} are small. In this 
regime we use the approximation 



T(W) ~- (,-.0) (A.5) 



r(x-Hy) y 

Then c(k) is approximated by 



cW-WIM-**)" 1 ( A - 6 ) 

i=i j>fc 

(where we have left out an irrelevant proportionality constant). As was mentioned previ- 
ously, the zonal spherical functions are not exactly known in this case, but we can use the 
large-x asymptotic formula 



Mrf ~ J~ 1/2 (x) E °( rk ) ei{rk,x) ( A - 7 ) 

Here the sum is over all the configurations obtained by the action of the elements of the 
Weyl group, i.e., the symmetry group of the restricted root system. The elements of the 
Weyl group are reflections and permutations. Since c(k) is invariant under reflections, only 
the even part J2 r c(rk) cos(rk, x) of J2 r c(rk) e l ( rk ' x ^ survives. Under permutations, 
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r=p 



c(rk) cos(rk,x) = l) CT ^c(/c) cos(k p Q)Xj) 



r=p 



z P (-iy {p) cos(k Pij)Xj : 

\.U=1 lij>k\^j K k) 



(A. 



Defining 



deti<ij< m [Cy = ^(-l) <j(p) cos(/Cp (i )X i ; 



deti< fc) K TO [A w ] = Y[{k 2 k - k 2 ) = det 1 < fe)/ < m /c 

k>l 



2(k-l) 



(A.9) 
(A.10) 



where detA is a Vandermonde determinant, we finally obtain from eqs. (jA.lj) . ()A.6|) . (I A. 71) 
and (TOl) 



P(x,s) oc Jj(x) / ] [ dk n k 2 n e~^ detAdetC 



n=l 



/m 
Y[ dk n deti^ijx, 
n=l 



.fc=i 



•'WEE - Z)(-ir W y dk k k 2 k e-~cos(k k x a )... J dk k ,k 2 k 7 l e~^T cos{k k ,x v ) 



a k=l k'=l 



(All) 



where we have used detldet-B = det[AB] and expanded the determinant according to 
detl = l) a ^ Ai a A2/3...A mu . We now use the known result for the integral 



" dA;fc 2n e- fe2fc2 cos(aA;) oc e-^# 2n (|- 



(A.12) 



where i^n denotes a Hermite polynomial of order 2n, to obtain 



m m 2 



y-,. _i 9 ) v... J2 ° ^~ u - 1 ' 

k=l k'=l 



...e 2s H- 



2m 



(A.13) 
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Setting 



(A.14) 



for brevity, the expression 



Y,{-iy {p) H 2 {x a )...H 2m {x v ) 



(A.15) 



is a determinant of Hermite polynomials of even order. The general form of H 2n (x) is 



H 2 „(x) = co(x 2n + Cl x 2n -' 2 + ... + c n ) 



2n-2 



(A.16) 



where c n ^ 0. To evaluate the determinant, we take out the irrelevant constants Co in front 
of the Hermite polynomials from the determinant so that we have monic polynomials in 
the latter. By adding an appropriate multiple of the first row to each of the other rows, we 
can eliminate the constants c n from each matrix element, except for the ones in the first 
row. Using that H 2 (x) = 4(x 2 — |) the determinant then takes the form 



deti< i)j < m [f/'2»(^)] oc 



z.2 1 



X 



x 2 p 2 (xi y 



r 2 — I 

x 2 2 



x 2 p 2 {x 2 ) 

£lP2(m-l)(£l) xlp 2 (m-i)(x 2 ) 



T 2 — I 

•^m 2 



X 2 m P 2 (x r 



-2 ~2 

c 1 ...x m 



P2(Xl) P 2 {X 2 ) 



xlpiixx) x\p 2 {x 2 ) 

x{p 2 (rn-l){xi) xlp 2 (rn-l){x 2 ) 
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P2{x m ) 



P2(m-l){Xl) P 2 (m-l)(x 2 ) ••• P 2 {m-l){x m ) 



+ 



X 2 m P 2 (Xm) 



x 2 m p 2 { m -i){x m ) 



(A.17) 



where p2 n denotes a polynomial of degree 2n and we have used general properties of deter- 
minants. Since we are interested in the insulating regime where 1 <C X\ <C £2 ••• <C x m , 
we can neglect the second determinant in the last line of the equation, whereas the first 
determinant can be re-written as a Vandermonde determinant by simply writing each line 
of the determinant as a linear combination of the other rows. Then we obtain 



P({xi}, s) oc sj J(x) J[ e 2 l s det!< itj < 
i=i 

m m m x a m m 

= [] sinh r?/2 (2x i ) J] smh^ 2 { Xj - x k )smhf 3/2 {x 3 + x k ) J[ sinh 9/2 x,e-^ JJ(xJ - x\) J] 4 

i=l j>k 1=1 P>q r=l 

(A.18) 

In the last line we have used ()5.6|) with mi, m , m s instead of the explicit root multiplicities. 
After using the identity for hyperbolic functions 

sinh(xj — Xjt)sinh(a;j + x k ) = sinh 2 ^- — smli 2 ^ (A. 19) 

()A.18|) becomes the expression reported in eq. (|6.6j) . 
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